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Abstract. A "pairing function" J associates a unique natural number z 
to any two natural numbers x,y such that for two "unpairing functions" 
K and L, the equalities K(J(x,y))=x, L(J(x,y))=y and J(K(z),L(z))=z 
hold. Using pairing functions on natural number representations of truth 
tables, we derive an encoding for Binary Decision Diagrams with the 
unique property that its boolean evaluation faithfully mimics its struc- 
tural conversion to a a natural number through recursive application 
of a matching pairing function. We then use this result to derive rank- 
ing and unranking functions for BDDs and reduced BDDs. The paper 
is organized as a self-contained literate Prolog program, available at 
http : //logic . csci . unt . edu/tarau/research/2008/pBDD . zip 
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1 Introduction 

This paper is an exploration with logic programming tools of ranking and un- 
ranking problems on Binary Decision Diagrams. The practical expressiveness of 
logic programming languages (in particular Prolog) are put at test in the pro- 
cess. The paper is part of a larger effort to cover in a declarative programming 
paradigm, arguably more elegantly, some fundamental combinatorial generation 
algorithms along the lines of pQ. However, our main focus is by no means "yet an- 
other implementation of BDDs in Prolog" . The paper is more about fundamental 
isomorphisms between logic functions and their natural number representations, 
in the tradition of [2], with the unusual twist that everything is expressed as a 
literate Prolog program, and therefore automatically testable by the reader. One 
could put such efforts under the generic umbrella of an emerging research field 
that we would like to call executable theoretical computer science. Nevertheless, 
we also hope that the more practically oriented reader will be able to benefit 
from this approach by being able to experiment with, and reuse our Prolog code 
in applications. 

The paper is organized as follows: Sections [2J and [3] overview efficient eval- 
uation of boolean formulae in Prolog using bitvectors represented as arbitrary 
length integers and Binary Decision Diagrams (BDDs). 



Section 2] discusses classic pairing and unpairing operations and introduces 
pairing/unpairing predicates acting directly on bitlists. 

Section [5] introduces a novel BDD encoding (based on our unpairing func- 
tions) and discusses the surprising equivalence between boolean evaluation of 
BDDs and the inverse of our encoding, the main result of the paper. 

Section [5] describes ranking and unranking functions for BDDs and reduced 
BDDs. 

Sections [7] and [5] discuss related work, future work and conclusions. 
The code in the paper, embedded in a literate programming LaTcX file, is 
entirely self contained and has been tested under SWI-Prolog. 

2 Parallel Evaluation of Boolean Functions with 
Bitvector Operations 

Evaluation of a boolean function can be performed one value at a time as in the 
predicate if _then_else/4 

if.then.elseCX.Y.Z.R) :- 
bit(X) ,bit(Y) ,bit(Z) , 
( X=1->Rf=Y 
; R=Z 
). 

bit(O) . 
bit(l) . 

resulting in a truth tabl^\ 

?- if_then_else(X,Y,Z,R) , write ([X.Y.Z] :R) ,nl,fail;nl. 

[0, 0, 0] :0 

[0, 0, 1] :1 

[0, 1, 0] :0 

[0, 1, 1] :1 

[1, 0, 0] :0 

[1, 0, 1] :0 

[1, 1, 0] :1 

[1, 1, 1]:1 

Clearly, this does not take advantage of the ability of modern hardware to per- 
form such operations one word a time - with the instant benefit of a speed-up 
proportional to the word size. An alternate representation, adapted from pQ uses 
integer encodings of 2" bits for each boolean variable Xq, . . . , X n _i. Bitvector 
operations evaluate all value combinations at once. 

1 One can see that if the number of variables is fixed, we can ignore the bitsrings 
in the brackets. Thus, the truth table can be identified with the natural number, 
represented in binary form by the last column. 



Proposition 1 Let x k be a variable for < k < n where n is the number of 
distinct variables in a boolean expression. Then column k in the matrix represen- 
tation of the inputs in the the truth table represents, as a bitstring, the natural 
number: 

x k = (2 2 " - l)/(2 2 """ + 1) (1) 

For instance, if n = 2, the formula computes xq = 3 = [0, 0, 1, 1] and x\ = 5 = 
[0,1,0,1]. 

The following predicates, working with arbitrary length bitstrings are used 
to evaluate variables Xk with k £ [0..n — 1] with formula[T]and map the constant 
boolean function 1 to the bitstring of length 2™, 111 . .1, representing 2 2 — 1 

•/. maps variable K in [0 . . NbOf Bits-1] to Xk 
var_to_bitstring_int (NbOf Bit s , K , Xk) : - 

all_ones_mask(Nb0f Bits, Mask) , 

var_to_bitstring_int (NbOf Bit s , Mask , K , Xk) . 

var_to_bitstring_int (NbOf Bit s , Mask , K , Xk) : - 
NK is NbOf Bits- (K+l) , 
D is (l«(l«3JK))+l, 
Xk is Mask//D. 

7. represents constant 1 as 11... 1 build with NbOf Bits bits 
all_ones_mask (NbOf Bits, Mask) : -Mask is (l«(l«Nb0f Bits) ) -1 . 

We have used in var_to_bitstring_int an adaptation of the efficient bitstring- 
integer encoding described in the Boolean Evaluation section of PQ . Intuitively, it 
is based on the idea that one can look at n variables as bitstring representations 
of the n columns of the truth table. 

Variables representing such bitstring-truth tables (seen as projection func- 
tions) can be combined with the usual bitwise integer operators, to obtain new 
bitstring truth tables, encoding all possible value combinations of their argu- 
ments. Note that the constant is represented as while the constant 1 is 
represented as 2 2 — 1, corresponding to a column in the truth table containing 
ones exclusively. 

3 Binary Decision Diagrams 

We have seen that Natural Numbers in [0..2 2 " — 1] can be used as represen- 
tations of truth tables defining n- variable boolean functions. A binary decision 
diagram (BDD) [3] is an ordered binary tree obtained from a boolean function, 
by assigning its variables, one at a time, to (left branch) and 1 (right branch). 
In virtually all practical applications BDDs are represented as DAGs after de- 
tecting shared nodes. We safely ignore this here as they represent the same logic 
function, which is all we care about at this point. Typically in the early litera- 
ture, the acronym ROBDD is used to denote reduced ordered BDDs. Because 



this optimization is now so prevalent, the term BDD is frequently use to refer 
to ROBDDs. Strictly speaking, BDD in this paper will stand for ordered BDD 
with reduction of identical branches but without node sharing. 

The construction deriving a BDD of a boolean function / is known as Shan- 
non expansion [I] , and is expressed as 

/(i) = (iA/[i*-0])V(iA/^l]) (2) 

where f[x<—a] is computed by uniformly substituting a for x in /. Note that 
by using the more familiar boolean if-the-else function Shannon expansion can 
also be expressed as: 

f{x) — if x then f[x *— 1] else f[x <— 0] (3) 

We represent a BDD in Prolog as a binary tree with constants and 1 as 
leaves, marked with the function symbol c/1. Internal if-then-else nodes marked 
with ite/3 are controlled by variables, ordered identically in each branch, as 
first arguments of ite/1. The two other arguments are subtrees representing 
the Then and Else branches. Note that, in practice, reduced, canonical DAG 
representations are used instead of binary tree representations. 

Alternatively, we observe that the Shannon expansion can be directly derived 
from a 2™ size truth table, using bitstring operations on encodings of its n vari- 
ables. Assuming that the first column of a truth table corresponds to variable 
x, x — and x — 1 mask out, respectively, the upper and lower half of the truth 
table. 

'/, splits a truth table of NV variables in 2 tables of NV-1 variables 
shannon_split (NV,X, Hi ,Lo) : - 

all_ones_mask(NV,M) , 

NV1 is NV-1, 

all_ones_mask (NV1 , LM) , 

HM is xor(M.LM) , 

Lo is A(LM,X), 

H is A(HM,X) , 

Hi is H»(1«NV1) . 

Note that the operation shannon_split can be reversed as follows: 

'/, fuses 2 truth tables of NV-1 variables into one of NV variables 
shannon_f use (NV , Hi , Lo , X) : - 

NV1 is NV-1, 

H is Hi«(l«NVl) , 

X is \/(H,Lo) . 

?- shannon_split (2, 7, X, Y) , shannon_fuse (2 , X,Y, Z) . 
X = 1, 
Y = 3, 
Z = 7. 

?- shannon_split (3, 42, X , Y) , shannon_fuse (3, X,Y, Z) . 



X = 2, 
Y = 10, 
Z = 42. 

Another way to look at these two operations (for a fixed value of NV), is 
as bijections associating a pair of natural numbers to a natural number, i.e. as 
pairing functions. 

4 Pairing and Unpairing Functions 

Definition 1 A pairing function is a bijection f : Nat x Nat — > Nat. An 
unpairing function is a bijection g : Nat — > Nat x Nat. 

Following Julia Robinson's notation [5], given a pairing function J, its left 
and right inverses K and L are such that 

J{K{z),L{z)) = z (4) 
K{J{x,y))=x (5) 

L(J(x,y))=y (6) 

We refer to 6 for a typical use in the foundations of mathematics and to 
[7] for an extensive study of various pairing functions and their computational 
properties. 

4.1 Cantor's Pairing Function 

Starting from Cantor's pairing function 

cantor_pair(Kl,K2,P) :-P is (((Kl-f*2)*(Kl-+K2+l))//2)lK2. 

bijections from Nat x Nat to Nat have been used for various proofs and con- 
structions of mathematical objects |5l6j . 

For X, Y £ {0, 1, 2, 3} the sequence of values of this pairing function is: 

?- f indalKR, (between(0,3,A) ,between(0,3,B) , cantor_pair (A,B,R) ) ,Rs) . 
Rs = [0, 2, 4, 6, 1, 5, 9, 13, 3, 11, 19, 27, 7, 23, 39, 55] 

Note however, that the inverse of Cantor's pairing function involves potentially 
expensive floating point operations that are also likely to loose precision for 
arbitrary length integers. 



4.2 The Pepis-Kalmar Pairing Function 

Another pairing function that can be implemented using only elementary integer 
operations is the following: 

/(as,y) = 2*(2y+l)-l (7) 

The predicates pepis_pair/3 and pepis_unpair/3 are derived from the function 
pepis_J and its left and right unpairing companions pepis_K and pepis_L that 
have been used, by Pepis, Kalmar and Robinson in some fundamental work on 
recursion theory, decidability and Hilbert's Tenth Problem in |8|9|10j : 

pepis_pair (X,Y,Z) :-pepis_J(X,Y,Z) . 

pepis_unpair (Z,X,Y) :-pepis_K(Z,X) ,pepis_L(Z,Y) . 

pepis_J(X,Y, Z):-Z is ( (1«X) * ( (Y«l)+1) ) -1 . 
pepis_K(Z, X):-Z1 is Z+l , two_s (Zl ,X) . 

pepis_L(Z, Y):-Z1 is Z+l ,no_two_s (Zl ,N) , Y is (N-l)»l. 

two_s(N,R) :-even(N) , ! ,H is N»l,two_s(H,T) ,R is T+l . 
two_s(_,0) . 

no_two_s(N,R) :-two_s(N,T) ,R is N // (1«T) . 

even(X) :- =:= A(1,X) . 
odd(X) :- 1 =:= A(1,X) . 

This pairing function is asymmetrically growing (faster growth on the first ar- 
gument). It works as follows: 

?- pepis_pair (1, 10,R) . 
R = 41. 

?- pepis_unpair ( 10 , 1 , R) • 
R = 3071. 

?- findall(R, (between (0,3, A) ,between(0,3,B) ,pepis_pair (A,B,R)) ,Rs) . 
Rs=[0, 2, 4, 6, 1, 5, 9, 13, 3, 11, 19, 27, 7, 23, 39, 55] 



4.3 Pairing/Unpairing operations acting directly on bitlists 

We will describe here pairing operations, that are expressed exclusively as bitlist 
transformations of bitmerge_unpair and its inverse bitmerge_pair, and are 
therefore likely to be easily hardware implementable. As we have found out 
recently, they turn out to be the same as the functions defined in Steven Pigeon's 
PhD thesis on Data Compression [IT], page 114). 

The predicate bitmerge_pair implements a bijection from Nat x Nat to Nat 
that works by splitting a number's big endian bitstring representation into odd 



and even bits, while its inverse to_pair blends the odd and even bits back to- 
gether. The helper predicates to_rbits and f romjrbits, given in the Appendix, 
convert to/from integers to bitlists. 

bitmerge_pair (X , Y , P) : - 
to_rbits(X,Xs) , 
to_rbits(Y,Ys) , 
bitmix(Xs,Ys,Ps) , ! , 
from_rbits(Ps,P) . 

bitmerge_unpair (P , X , Y) : - 
to_rbits(P,Ps) , 
bitmix(Xs,Ys,Ps) , ! , 
f rom_rbits (Xs , X) , 
from_rbits(Ys,Y) . 

bitmix([X|Xs] ,Ys, [X|Ms]) :-! ,bitmix (Ys ,Xs ,Ms) . 
bitmix([] , [X|Xs] , [0|Ms]) :-! ,bitmix( [X | Xs] , [] ,Ms) . 
bitmix([] ,[],[]). 

The transformation of the bitlists, done by the bidirectional predicate bitmerge 
is shown in the following example with bitstrings aligned: 

?- bitmerge_unpair (2008 ,X ,Y) ,bitmerge_pair (X,Y,Z) . 
X = 60, 
Y = 26, 
Z = 2008 

•/„ 2008: [0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1] 
I 60: [ 0, 1, 1, 1, 1] 

I 26: [ 0, 1, 0, 1, 1 ] 

Note that we represent numbers with bits in reverse order (least significant on 
the left). Like in the case of Cantor's pairing function, we can see similar growth 
in both arguments: 

?- between (0 , 15 , N) , bitmerge_unpair (N , A , B) , 

write(N: (A,B) ) ,write(' ') ,fail;nl. 
0: (0, 0) 1: (1, 0) 2: (0, 1) 3: (1, 1) 
4: (2, 0) 5: (3, 0) 6: (2, 1) 7: (3, 1) 
8: (0, 2) 9: (1, 2) 10: (0, 3) 11: (1, 3) 
12: (2, 2) 13: (3, 2) 14: (2, 3) 15: (3, 3) 

?- between(0,3,A) , between (0, 3, B) ,bitmerge_pair (A,B,N) , 

write(N: (A,B) ) , write ( ' ') ,fail;nl. 
0: (0, 0) 2: (0, 1) 8: (0, 2) 10: (0, 3) 
1: (1, 0) 3: (1, 1) 9: (1, 2) 11: (1, 3) 
4: (2, 0) 6: (2, 1) 12: (2, 2) 14: (2, 3) 
5: (3, 0) 7: (3, 1) 13: (3, 2) 15: (3, 3) 

It is also convenient sometimes to see pairing/unpairing as one-to-one functions 
from/to the underlying language's ordered pairs, i.e. X-Y in Prolog : 



bitmerge_pair (X-Y,Z) :-bitmerge_pair (X,Y,Z) . 
bitmerge_unpair (Z,X-Y) : -bitmerge_unpair (Z, X, Y) . 



5 Encodings of Binary Decision Diagrams 

We will build a BDD by applying bitmerge_unpair recursively to a Natural 
Number TT, seen as an TV-variable 2 N bit truth table. This results in a complete 
binary tree of depth N. As we will show later, this binary tree represents a BDD 
that returns TT when evaluated applying its boolean operations. 

'/. NV=number of varibles, TT=a truth table, BDD the result 
plain_bdd (NV , TT , bdd (NV , BDD) ) : - 

Max is (1«(1«NV)) , 

TT<Max , 

isplit(NV,TT, BDD) . 

"/, recurses to depth NV, splitting TT into pairs 
isplit(0,TT,c(TT)) . 
isplit(NV,TT,R) :-NVX), 
NV1 is NV-1, 

bitmerge_unpair (TT,Hi,Lo) , 
isplit(NVl,Hi,H) , 
isplit(NVl,Lo,L) , 
ite(NVl,H,L)=R. 

The following examples show the results returned by plain_bdd for all 2 2> " truth 
tables associated to k variables, with k = 2. 

?- between (0, 15, TT) ,plain_bdd(2,TT,BDD) , write (TT: BDD) ,nl,fail;nl 
0:bdd(2, ite(l, ite(0, c(0), c(0)), ite(0, c(0), c(0)))) 
l:bdd(2, ite(l, ite(0, c(l), c(0)), ite(0, c(0), c(0)))) 
2:bdd(2, ite(l, ite(0, c(0), c(0)), ite(0, c(l), c(0)))) 

13:bdd(2, ite(l, ite(0, c(l), c(l)), ite(0, c(0), c(l)))) 
14:bdd(2, ite(l, ite(0, c(0), c(l)), ite(0, c(l), c(l)))) 
15:bdd(2, ite(l, ite(0, c(l), c(l)), ite(0, c(l), c(l)))) 



5.1 Reducing the BDDs 

The predicate bdd_reduce reduces a BDD by trimming identical left and right 
subtrees, and the predicate bdd associates this reduced form to N G Nat. 

bdd_reduce(BDD,bdd(NV,R)) :-nonvar(BDD) ,BDD=bdd(NV,X) ,bdd_reducel (X,R) . 
bdd_reducel(c(TT) ,c(TT)) . 

bdd_reducel (ite(_,A,B) ,R) : -A=B,bdd_reducel (A,R) . 
bdd_reducel(ite(X,A,B) ,ite(X,RA,RB)) :-A\=B, 



bdd_reducel (A,RA) ,bdd_reducel (B,RB) . 

bdd(NV,TT, ReducedBDD) : - 
plain.bddCNV.TT, BDD) , 
bdd_reduce (BDD , ReducedBDD) . 

Note that we omit here the reduction step consisting in sharing common subtrees, 
as it is obtained easily by replacing trees with DAGs. The process is facilitated 
by the fact that our unique encoding provides a perfect hashing key for each 
subtree. The following examples show the results returned by bdd for NV=2. 

?- between (0, 15 ,TT) ,bdd(2,TT,BDD) , write (TT: BDD) ,nl,fail;nl 
0:bdd(2, c(0)) 

l:bdd(2, ite(l, ite(0, c(l), c(0)), c(0))) 
2:bdd(2, ite(l, c(0), ite(0, c(l), c(0)))) 
3:bdd(2, ite(0, c(l), c(0))) 

13:bdd(2, ite(l, c(l), ite(0, c(0), c(l)))) 
14:bdd(2, ite(l, ite(0, c(0), c(D), c(l))) 
15:bdd(2, c(l)) 



5.2 From BDDs to Natural Numbers 

One can "evaluate back" the binary tree representing the BDD, by using the 
pairing function bitmerge_pair. The inverse of plain_bdd is implemented as 
follows: 

plain_inverse_bdd(bdd(_,X) ,TT) : -plain_inverse_bddl (X,TT) . 

plain_inverse_bddl (c(TT) , TT) . 
plain_inverse_bddl (ite(_,L,R) , TT) : - 

plain_inverse_bddl (L,X) , 

plain_inverse_bddl (R,Y) , 

bitmerge_pair (X,Y,TT) . 

?- plain_bdd(3,42, BDD) ,plain_inverse_bdd(BDD,N) . 
BDD = bdd (3, 

ite(2, 

ite(l, 

ite(0, c(0), c(0)), 
ite(0, c(0), c(0))), 
ite(l, 

ite(0, c(l), c(l)), 
ite(0, c(l), c(0))))), 

N = 42 

Note however that plain_inverse_bdd/2 does not act as an inverse of bdd/3, 
given that the structure of the BDD tree is changed by reduction. 



5.3 Boolean Evaluation of BDDs 



This raises the obvious question: how can we recover the original truth table from 
a reduced BDD? The obvious answer is: by evaluating it as a boolean function! 
The predicate ev/2 describes the BDD evaluator: 

ev(bdd(NV,B) ,TT) : - 
all_ones_mask(NV,M) , 
eval_with_mask (NV,M,B,TT) . 

evc(0,_,0) . 
evc(l,M,M) . 

eval_with_mask(_,M,c(X) ,R) :-evc(X,M,R) . 
eval_with_mask (NV , M , ite (X , T , E) , R) : - 

eval_with_mask (NV , M , T , A) , 

eval_with_mask (NV ,M,E,B) , 

var_to_bitstring_int (NV.M.X.V) , 

ite(V,A,B,R) . 

The predicate ite/4 used in eval_with_mask implements the boolean function 
if X then T else E using arbitrary length bitvector operations: 

ite(X,T,E, R):-R is xor(/\(X,xor(T,E)) ,E) . 

Note that this equivalent formula for ite is slightly more efficient than the 
obvious one with A and V as it requires only 3 boolean operations. We will 
use ite/4 as the basic building block for implementing a boolean evaluator for 
BDDs. 

5.4 The Equivalence 

A surprising result is that boolean evaluation and structural transformation with 
repeated application of pairing produce the same result, i.e. the predicate ev/2 
also acts as an inverse of bdd/2 and plain_bdd/2. 

,4s the following example shows, boolean evaluation ev/2 faithfully emulates 
plain_inverse_bdd/2, on both plain and reduced BDDs. 

?- plain_bdd(3,42,BDD) ,ev(BDD,N) . 
BDD = bdd(3, 

ite(2, 

ite(l, 

ite(0, c(0), c(0)), 
ite(0, c(0), c(0))), 
ite(l, 

ite(0, c(l), c(l)), 
ite(0, c(l), c(0))))), 

N = 42 

?- bdd(3,42,BDD) ,ev(BDD,N) . 
BDD = bdd(3, 



ite(2, 
c(0), 
iteCl, 

c(l), 

ite(0, c(l), c(0))))), 

N = 42 

The main result of this subsection can now be summarized as follows: 

Proposition 2 Let B be the complete binary tree of depth N , obtained by re- 
cursive applications o/bitmerge_unpair on a truth table T, as described by the 
predicate plain_bdd(N,T,B) . 

Then for any N and any T , when B is interpreted as an (unreduced) BDD, 
the result V of its boolean evaluation using the predicate ev(N, B, V) and the 
result R obtained by applying plain JnverseJbdd(N, B, R) are both identical to 
T. Moreover, the operation ev(N, B,V) reverses the effects of both plain_bdd 
and bdd with an identical result. 

Proof: The predicate plain_bdd builds a binary tree by splitting the bitstring 
tt E [Q..2 N — 1] up to depth N. Observe that this corresponds to the Shannon 
expansion [4J of the formula associated to the truth table, using variable order 
[n — 1, 0]. Observe that the effect of bitstring_unpair is the same as 

— the effect of var_to_bitstring_int(N,M, (N-l) ,R) acting as a mask select- 
ing the left branch 

— and the effect of its complement, acting as a mask selecting the right branch. 

Given that 2 N is the double of 2 JV ~ 1 , the same invariant holds at each step, as 
the bitstring length of the truth table reduces to half. On the other hand, it 
is clear that ev reverses the action of both plain_bdd and bdd as BDDs and 
reduced BDDs represent the same boolean function [3]. 

This result can be seen as a yet another intriguing isomorphism between 
boolean, arithmetic and symbolic computations. 

6 Ranking and Unranking of BDDs 

One more step is needed to extend the mapping between BDDs with N variables 
to a bijective mapping from/to Nat: we will have to "shift toward infinity" the 
starting point of each new block of BDDs in Nat as BDDs of larger and larger 
sizes are enumerated. 

First, we need to know by how much - so we compute the sum of the counts 
of boolean functions with up to N variables. 

bsum(0,0) . 

bsum(N,S) :-NX),Nl is N-l ,bsuml (Nl ,S) . 
bsuml(0,2) . 

bsuml(N,S) :-NX),Nl is N-l ,bsuml (Nl , SI) , S is S1+(1«(1«N) ) . 



The stream of all such sums can now be generated as usual: 
bsum(S) :-nat(N) ,bsum(N,S) . 

nat (0) . 

nat(N) :-nat(Nl) ,N is Nl+1 . 

What we are really interested in, is decomposing N into the distance to the last 
bsum smaller than N, N_M and the index of that generates the sum, K. 

to_bsum(N, X,N_M) :- 

nat(X) ,bsum(X,S) ,S>N, ! , 
K is X-l, 
bsum(K,M), 
N_M is N-M. 

Unranking of an arbitrary BDD is now easy - the index K determines the number 
of variables and N_M determines the rank. Together they select the right BDD 
with plain_bdd and bdd/3. 

nat2plain_bdd(N,BDD) :-to_bsum(N, K,N_M) ,plain_bdd(K,N_M,BDD) . 
nat2bdd(N,BDD) :-to_bsum(N, K,N_M) ,bdd(K,N_M,BDD) . 

Ranking of a BDD is even easier: we first compute its NumberOf Vars and its rank 
Nth, then we shift the rank by the bsums up to NumberOf Vars, enumerating the 
ranks previously assigned. 

plain_bdd2nat (bdd (Number Of Vars ,BDD) ,N) :- 
B=bdd (Number Of Vars ,BDD) , 
plain_inverse_bdd(B,Nth) , 
K is NumberOf Vars-1 , 
bsum(K,S) ,N is SH4Jth. 

bdd2nat (bdd (NumberOf Vars , BDD) ,N) :- 
B=bdd (Number Of Vars ,BDD) , 
ev(B.Nth) , 

K is NumberOf Vars-1 , 
bsum(K,S),N is S+Mth. 

As the following example shows, nat2plain_bdd/2 and plain_bdd2nat/2 im- 
plement inverse functions. 

?- nat2plain_bdd(42,BDD) ,plain_bdd2nat (BDD.N) . 
BDD = bdd (4, 

ite(3, 

ite(2, 

ite(l, 

ite(0, c(0), c(0)), 
ite(0, c(l), c(0))), 
ite(l, 

ite(0, c(l), c(0)), 
ite(0, c(0), c(0)))), 



ite(2, 

iteCl, 

ite(0, c(0), c(0)), 
ite(0, c(0), c(0))), 
ite(l, ite(0, c(0), c(0)), 

ite(0, c(O), c(O)))))), 

N = 42 

The same applies to nat2bdd/2 and its inverse bdd2nat/2. 

?- nat2bdd(42,BDD) , bdd2nat (BDD , N) . 
BDD = bdd(4, 

ite(3, 

ite(2, 

ite(l, c(0), 

ite(0, c(l), c(0))), 
ite(l, 

ite(0, c(l),c(0)), c(0))), 
c(0))), 

N = 42 

We can now generate infinite streams of BDDs as follows: 
plain_bdd(BDD) :-nat(N) ,nat2plain_bdd(N,BDD) . 

bdd(BDD) :-nat(N) ,nat2bdd(N,BDD) . 

?- plain_bdd(BDD) . 

BDD = bdd(l, ite(0, c(0), c(0))) ; 

BDD = bdd(l, ite(0, c(l), c(0))) ; 

BDD = bdd(2, iteCl, ite(0, c(0), c(0)), ite(0, c(0), c(0)))) 

BDD = bdd(2, iteCl, ite(0, c(l), c(0)), ite(0, c(0), c(0)))) 

?- bdd(BDD) . 

BDD = bdd(l, c(0)) ; 

BDD = bdd(l, ite(0, c(l), c(0))) ; 

BDD = bdd(2, c(0)) ; 

BDD = bdd(2, iteCl, ite(0, c(l), c(0)), c(0))) ; 

BDD = bdd(2, iteCl, c(0), ite(0, c(l), c(0)))) ; 

BDD = bdd(2, ite(0, c(l), c(0))) ; 



7 Related work 

Pairing functions have been used in work on decision problems as early as |8I9I5] . 
Ranking functions can be traced back to Godel numberings [2112] associated to 
formulae. Together with their inverse unranking functions they are also used in 
combinatorial generation algorithms [1311] . Binary Decision Diagrams are the 
dominant boolean function representation in the held of circuit design automa- 
tion [Tlj ■ BDDs have been used in a Genetic Programming context [15116] as a 



representation of evolving individuals subject to crossovers and mutations ex- 
pressed as structural transformations and recently in a machine learning context 
for compressing probabilistic Prolog programs |17j representing candidate the- 
ories. Other interesting uses of BDDs in a logic and constraint programming 
context are related to representations of finite domains. In [18] an algorithm for 
finding minimal reasons for inferences is given. 

8 Conclusion and Future Work 

The surprising connection of pairing/unpairing functions and BDDs, is the in- 
direct result of implementation work on a number of practical applications. Our 
initial interest has been triggered by applications of the encodings to combina- 
tional circuit synthesis in a logic programming framework 19 20 . We have found 
them also interesting as uniform blocks for Genetic Programming applications 
of Logic Programming. In a Genetic Programming context [21], the bijections 
between bitvectors/natural numbers on one side, and trees/graphs representing 
BDDs on the other side, suggest exploring the mapping and its action on vari- 
ous transformations as a phenotype-genotype connection. Given the connection 
between BDDs to boolean and finite domain constraint solvers it would be inter- 
esting to explore in that context, efficient succinct data representations derived 
from our BDD encodings. 
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Appendix 

To make the code in the paper fully self contained, we list here some auxiliary 
functions. 

'/, converts an int to a list of bits, least significant first 
to_rbits(0, []) . 

to_rbits(N, [B | Bs] ) :-N>0,B is N mod 2, Nl is N//2, 
to_rbits(Nl,Bs) . 

'/, converts a list of bits (least significant first) into an int 
from_rbits(Rs,N) :-nonvar(Rs) , f rom_rbits (Rs , 0, ,N) . 

from_rbits([] ,_,N,N) . 

from_rbits([X|Xs] ,E,N1,N3) :-NewE is E+1.N2 is X«E-rNl, 
from_rbits(Xs,NewE,N2,N3) . 



